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Abstract 

The ground state entanglement of the two- mode Bose-Einstein condensate is investigated through 
a quantum phase transition approach. The entanglement measure is taken as the order parameter 
and this is a non-local order parameter, which is different from the conventional order parameter 
of the Mott insulator-superfluid phase transitions. For this non-local order parameter, scaling 
behavior corresponding to a continuous phase transition is obtained and a power-law divergence 
near the critical region follows it. This scaling behavior of quantum entanglement is analyzed 
by the finite-size scaling and the critical exponents are obtained as v = 1.01 and 7 = 0.86. A 
close connection between quantum fluctuations and the phase transition of entanglement is also 
obtained. 
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I. INTRODUCTION 

Quantum entanglement is a key feature of quantum information theory and it is one of the 
most active research areas in recent years [I], especially in the area of its combination with 
condensed matter systems [2]. Beyond its generation and application, one of the essential 
questions is how to understand the process that happens in a physical system when it transits 
from non-entangled states to entangled states. One approach to study this phenomenon is 
to treat it as a quantum phase transition, where the order parameter is the entanglement 
measure. This approach combines the theory of critical phenomenon with the theory of 
quantum information. For spin lattice models, many results have been obtained. The 
concurrence [31 H] is often used as the entanglement measure in spin models. Entanglement 
exhibits scaling behavior near the critical region and the critical behavior is shown to be 
logarithm [5HH]. The phase transition is second-order for the ferromagnetic case and is 
first-order for the antiferromagnetic case [PJ. 

While the critical behavior of quantum entanglement in spin models have widely been 
studied, there are very few studies on that in boson systems. It is thus of interest to 
investigate the critical behavior of quantum entanglement in boson systems. 

One extensively studied boson system in quantum entanglement is the two-mode Bose- 
Einstein condensates coupled via Josephson tunneling [TU]. It is described by the Hamilto- 
nian [TT] 



where ai,a 2 are the annihilation operators for the two modes (1 and 2) respectively, and 
N% = a\a\,N2 = a\a2 are the corresponding number operators. The parameter K provides 
the atom-atom interaction, Ap is the difference in the chemical potential between the two 
modes and Ej is the coupling for tunneling. This Hamiltonian describes both the double-well 
Bose-Einstein condensate and the two-level Bose-Einstein condensate in a single potential. 
For the first case, the tunneling between the two wells must be small to use this Hamiltonian, 
while for the second case, there is no such restriction. We will show in this paper that the 
phase transition occurs at very small couplings, so the quantum phase transition approach 
can describe both cases. The entanglement production in this system has been extensively 
studied [10J [T2T[T5] . The von Neumann entropy [16] E(p) is the usually used entanglement 
measure, where p is the density matrix of the system, and for a system size of iV particles, 
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the maximum entropy is E max = log 2 (iV +1). 

This Hamiltonian ([I]) is in fact a two-site version of the Bose- Hubbard model p2] . When 
varying the ratio between the interaction term the coupling term through a critical value, a 
quantum phase transition occurs in the Bose-Hubbard model, which is the Mott insulator to 
the superfluid transition [IB]. This phase transition is driven by quantum fluctuations and 
the order parameter is the conventional wave function. In the Mott insulator phase, atoms 
are localized in lattice sites, while in the superfluid phase, atoms spread out over the whole 
system. Although the insulator-superfluid phase transition is studied extensively [T9H23] 
both in theory and in experiment, it is interesting to investigate what would happen to the 
Bose-Hubbard model when taking a non-local order parameter, rather than the conventional 
order parameter. 

In this paper, we present such a study for the simplest two-site Bose-Hubbard system, 
the two-mode Bose-Einstein condensate. The entanglement measure, the von Neumann 
entropy, is taken as the non-local order parameter. We show that there is a critical point 
and entanglement exhibits scaling behavior near the critical point, which can be analyzed 
using the theory of critical phenomena. We identify this as a continuous phase transition. 
This phase transition is different from the insulator-superfluid phase transition, because it 
is obtained for a non-local order parameter, rather than the conventional order parameter. 
The critical behavior of quantum entanglement is shown to be power- law divergent, which 
is different from the logarithm divergent of spin lattice models. Our work may help the 
combination of methods in critical phenomena and quantum information theory for the 
boson systems, especially for the Bose-Hubbard model. Further extension of this work is to 
investigate quantum phase transitions in the Bose-Hubbard model of dozens of lattice site, 
where a new entanglement measure is also needed to be put up. 



II. CONTINUOUS PHASE TRANSITION 

In this paper, we only consider the case K > 0, corresponding to a repulsive interaction 
between atoms. The total particle number is conserved and we set A/i = 0. Using the 
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angular momentum operators 

Jm = \{N 2 -N 1 ) 
Jx = ^(a\ a 2 + 4 a i) 
Jy = ^( a i a 2 - 4 a i) 
and neglecting constant terms, the Hamiltonian (jlj is rewritten as 

H = X J 2 z -nj x , (2) 

where x = K/2 and Q = Ej. As we are only interested in the ratio between the two 
competing energy terms, it is convenient to introduce the dimensionless parameter in 
the calculation, so the Hamiltonian can be reduced to 

n = j 2 z - nj x , (3) 

where we have redefined Q using the dimensionless parameter, i.e., fi/x ~~ ► ^- This di- 
mensionless coupling parameter can be viewed as an 'external field' by analogy with Ising 
models. We also define the dimensionless entropy E(p)/E max — > E(p) to make it easier to 
compare the results of different system sizes. We use numerical diagonalization to calculate 
[2H EE] the ground state entanglement and its susceptibility with respect to the external 
field a 

We first calculate the susceptibility dE J^ with respect to the coupling VL for various system 
sizes, which is shown in Fig. [TJ We see that there is a critical point fl m for each system 
size, where the susceptibility reaches its critical value dE J^ ■ The critical susceptibility 
dE J^ increases with the system size and would be divergent for an infinite system size 
that corresponds to the thermodynamic limit, which implies that this is a continuous phase 
transition where there is no discontinuity in the order parameter, as depicted by the inset 
for the system of N = 2700 particles. This will be verified further in section |IV| 

From Fig. [TJ the critical point Q m lies in the small coupling regime, which means the 
phase transition occurs shortly after the external field is switched on. We could easily figure 
that the critical value is Q c = for an infinite system size of the thermodynamic limit. When 
Q = 0, the Neumann entropy is zero and there is no two-mode entanglement in the system; 
When Q > 0, the Neumann entropy gets a finite value and entanglement is generated in 
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FIG. 1. (color online) The susceptibility d ^f l of the quantum entanglement with respect to 
the external field Q for system sizes of N = 240, 400, 700, 1000, 2100, 2700. A continuous phase 
transition occurs as the susceptibility diverges with the system size. The critical point TO where 



the susceptibility attains its maximum 



dE{p) 
dtt i 



lies in the Fock regime and this critical susceptibility 



diverges with the system size. The inset depicts the change of the order parameter-the ground state 
entanglement for the system of N = 2700, which increases continuously from zero. We choose even 
particle numbers because for odd particle numbers there is a degeneracy of the ground state when 

n = o. 



the system. That means the system transits from non-entangled states to entangled states, 
two essentially different states, once Q is switched on from 0, so the critical value is just 0. 



This will be verified further in section |III[ where we numerically fit the critical point and 
the critical susceptibility for various system sizes. The critical point Q m is well fitted to iV 
by choosing Q c = 0. 



III. POWER-LAW DIVERGENCE 



The well-behaved relationship between the critical point and the system size in Fig. [2] is 
not just a coincidence. Actually it manifests the scaling behavior of quantum entanglement 
for this quantum system, which is typical in critical phenomenon. From Fig. |2j we obtain 
the scaling relationship 

n m = 0.319225A- - 989062 (4) 



for the critical point and the scaling relationship 

0.393037A - 846662 . (5) 



dE(p) ri .- )(V -, (V) _ v .o.s | ( ,(,i,2 



for the critical susceptibility. The scaling behavior of the susceptibility is power-law diver- 
gent, in contrast to the logarithm divergence of spin lattice systems [5]. 

This power-law divergence of the susceptibility can be understood in the thermodynamic 
limit using a simple analysis. The basic idea is to truncate the Fock space of the system to 
just three basis states and use them to approximate the state of the system. The validity 
of this approximation lies in the fact that the critical point is Sl c = and the derealization 
process is very weak near this critical point, which means that the transitions between 
different basis states of the original Fock space are very weak and we can use the three 
most important basis states for approximation. This is verified at the end of the calculation 
in Eq. ([6]), where a power-law behavior of the susceptibility is obtained and the divergent 
exponent does not differ much from that of the numerical simulation. 

We choose the Fock space basis \Nx,N 2 > for the system, where Ni is the particle number 
in the first site and N% is the particle number on the second site. When Q = 0, the ground 
state is \N/2, N/2 > with energy E = 0, that is, the system is in a self-trapping state without 
particle tunneling between the two sites. As Q increases, the particles begin tunneling 
between the two sites and this derealization process connects different basis states, so the 
system is described by Yln=o c n \n, N—n >. The critical value is fl c = and the derealization 
process is very weak near this region, so we can truncate the Fock space of the system to 
just three basis states | N/2, N/2 >, | N/2 - 1, N/2 + 1 > and | N/2 + 1, N/2 - 1 >, then the 
state of the system is \if> >= c \N/2,N/2 > +c x \N/2- 1, N/2 + 1 > +c 2 | N/2 + 1, N/2 - 1 >, 
where we assume the coefficient q to be real numbers for simplicity. As the probabilities of 
tunneling between the two sites are equal, the coefficients C\ and c 2 are equal. Combining 



with the normalization condition c\ + c\ + c 2 = 1, we get the relationship Co = a/1 — 2cf. 
We next calculate the approximate ground state to determine the value of the coefficients, 



ttN N N -^fJl-2c\ + c Xn N N ,N N 

ww > = -— «{i T , Y > + 4V _^ Ci ' n T - 1, T + 1 > +i T + 1, T - 1 >» 

= BW >, 

where the approximation y / N/2(N/2 + 1) ~ iV/2 is taken. The critical point fi c = 
determines that C\ is a small number. From \ip >— Cq[\N/2, N/2 > +c\/cq{\N/2 — 1, N/2 + 
f > +\N/2 + 1, N/2 - 1 >)], we obtain 

2 c ' 

which is approximately zero and is the ground state energy near Q c = 0, and 
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which gives the value 



i- 1 
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There are two values of c\ and what we choose is the smaller one. Substituting the values 
of the coefficients into the von Neumann entropy 

E{p) = -cl log 2 eg - Ci log 2 Ci - c 2 2 log 2 c 2 2 

and taking its derivative with respect to Q gives 



dE(p) ttN 
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j N" M7 (6) 

where the relationship flN ~ jV - 01 from Eq. Q in the thermodynamics limit is used. Thus 
we briefly illustrate the power-law divergence of the susceptiblity in the thermodynamic 
limit. 

The divergent exponent obtained in the analytic calculation is 0.97 and it is different from 
the value 0.85 of the numerical simulation in Eq. ([5]). This difference may be accounted by 
the finite-size effects and the truncation errors. First, the analytic calculation manifests the 
thermodynamic limit, where there is no finite-size effect. The numerical result, however, 
is influenced by the finite-size effects, so this may be one of the reasons for the difference 
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FIG. 2. (color online) The scaling behavior of the quantum phase transition with the critical 
value O c = 0. The critical point O m approaches asymptotically by SI 

^-0.989062 and the 

susceptibility diverges in a power-law behavior captured by dE J^ ~ jy - 846662 , which is different 
from the logarithm divergence of spin lattice models. 



between the divergent exponents. Second, we adopt approximation in the analytic calcula- 
tion by truncating the Fock space of the system to just three basis states. The numerical 
simulation, however, includes the full Fock space. The neglected basis states would certainly 
contribute to the result, although their amplitudes are small near the critical point. So the 
difference between the divergent exponent is also influenced by the truncation errors. 

IV. FINITE-SIZE SCALING 

A key feature of the critical phenomenon is the finite-size scaling. Phase transitions 
only occur at the thermodynamic limit, while numerical simulations can only deal with 
finite system sizes. To extract information from results obtained from the finite system, 
the finite-size scaling is required, where the effect of finite system sizes are eliminated by 
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collecting all data of various system sizes onto a single curve and the critical exponent can 
be deduced in this process. In the phase transition of thermal order parameters, e. g. , the 
magnetization, the critical exponent v of the correlation length satisfies |T — T c \ ~ N~ l l u . 
By analogy, we obtain v = 1/0.989062 ~ 1.01 from Eq. [4| which is the critical exponent for 
the quantum phase transition of quantum entanglement. This critical exponent gives the 
reduced coordinate N u (Cl — Cl m ) for all the finite system sizes. From Eq. [5j the susceptibility 
is reduced to jv~ 0,85 ( d ^^ — dE J^ )■ If the quantum entanglement of this model manifests 
quantum phase transitions, then all data of various system sizes could be collected onto a 
single cure using the above reduced coordinates. This is indeed the case as exhibited in Fig. 
[3j Again this relationship is not just a coincidence. It illustrates that quantum entanglement 
of this model indeed belongs to critical phenomenon. Resorting to the the phase transition 
of the Magnetization, where the susceptibility \ °f the Magetization is reduced to N^^^x-, 
we obtain the critical exponent 7 = 0.85u ~ 0.86 in this model. 

V. QUANTUM FLUCTUATIONS 

The Mott insulator-superfluid transition is driven by quantum fluctuations, which is com- 
mon for quantum phase transitions. Here we show that a close connection also exists between 
quantum fluctuations and the phase transition of entanglement . In the dynamical regime of 
entanglement production, the system is required to undergo a derealization process, where 
large quantum fluctuation exists, to generate entanglement. So the quantum phase tran- 
sition of entanglement should be closely related to quantum fluctuations. In the angular 
momentum representation \j,j z >, where j = N/2 and j z = —N/2, —N/2 + 1, . . . , N/2, 
the quantum fluctuation is (AJ Z ) 2 = (J 2 ) — (J z ) 2 . We show that (AJ Z ) 2 and E(p) have a 
similiar behavior, which indicates their close connection with each other. We plot (AJ Z ) 2 , 
E(p) and their derivatives with respect to in their reduced value in Fig. |4j We see that 
both the quantum fluctuation and the quantum entanglement grows with the external field 
Q, and their growth corresponds to each other, which can be seen from their derivatives. As 
quantum entanglement is a non-classical correlation, it is consistent that its quantum phase 
transition is closely related to quantum fluctuations. 

There is a 'delay' between the derivative of E(p) and that of (A J z ) 2 , where the derivative 
of E(p) reaches its maximum value earlier than the derivative of (AJ 2 ) 2 . This is due to the 
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FIG. 3. (color online) The finite-size scaling for the quantum phase transition of quantum en- 
tanglement. After the susceptibility is reduced by the system size to jV~°- 85 ( rf ^'^ — ), it 
becomes a function of ./V 1,01 (J) — Q m )- Data from a broad range of system sizes are collected on 
this single curve. The critical exponent obtained is v = 1.01 and 7 = 0.86. 



finite-size effects. We are not working in the thermodynamic limit, so the derivatives between 
the quantum fluctuation and the quantum entanglement are not in complete correspondance. 
This is further confirmed by Table |TJ where the 'delay' AQ between the maximum points of 
the derivatives are calculated for various system sizes. We see that the 'delay' between them 
is decreasing as the system size growes, so we can figure that in the thermodynamic limit, 
the growth behavior of the entanglement and that of the fluctuation will approximately 
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FIG. 4. (color online) Comparison between (AJ 2 ) 2 , E(p) and their derivatives with respect to fL 
Their values are reduced to 1 by their maximum values. The increase of the quantum fluctuation 
corresponds to the increase of the order parameter, which indicates its connection with the phase 
transition of entanglement. There is a small 'delay' between the derivative of the fluctuation and 
the susceptibility, with the susceptibility obtaining maximum value first. This 'delay' comes from 
the finite-size effects. 



correspond to each other. 



VI. SUMMARY 

In summary, we have studied the entanglement of a boson system from the quantum 
phase transition approach. It is shown that in this system there is a continuous phase tran- 
sition for the non-local order parameter and entanglement exhibits scaling behavior near 
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TABLE I. The 'delay' AQ between the maximum points of the derivatives of E(p) and (AJ 2 ) 2 for 
various system sizes. As the system size increases, the 'delay' decreases, which means the growth 
behavior of the entanglement and that of the fluctuation are more closely related. This suggests 
that in the thermodynamic limit, the two growth behaviors will correspond to each other. 
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the critical point, with the critical exponent calculated to be u = 1.01 and 7 = 0.86. The 
critical behavior under discussion is different from that of the spin lattice models, because 
a power-law divergence is obtained for the boson system while it is logarithm divergence 
for the spin models. A further study of this phenomenon may consist of putting up an 
entanglement measure for Boson systems of more lattice sites, investigating the quantum 
phase transition of the Bose-Hubbard model of more lattice sites and obtaining its univer- 
sality class. The renormalization group method specially for taking into account the effect 
of quantum entanglement [26H28] may be used in that case. 
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